Domain wall structure in magnetic bilayers with perpendicular 
anisotropy 
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PACS 75 . 60 . Ch - Domain walls and domain structure 

PACS 75.70.-i - Magnetic properties of thin films, surfaces, and interfaces 

PACS 68 . 37 . Ef - Scanning tunneling microscopy (including chemistry induced with STM) 

Abstract. - We study the magnetic domain wall structure in magnetic bilayers (two ultrathin 
ferromagnetic layers separated by a non magnetic spacer) with perpendicular magnetization. Com- 
bining magnetic force and ballistic electron emission microscopies, we are able to reveal the details 
of the magnetic structure of the wall with a high spatial accuracy. In these layers, we show that 
the classical Bloch wall observed in single layers transforms into superposed Neel walls due to the 
magnetic coupling between the ferromagnetic layers. Quantitative agreement with micromagnetic 
calculations is achieved. 



Introduction. — Magnetic domain walls (DW) are 
attracting a sustained interest both from a fundamental 
perspective and for applications, especially in the context 
of innovative magnetic memories pji2j . The aim is to con- 
trol the DW motion using magnetic fields or spin polar- 
ized currents. The dynamics of these walls is intrinsically 
linked to their micromagnetic structures, which impact 
their mobility and the so-called Walker velocity (the crit- 
ical field or current above which the DW magnetization 
precesses when moving) [3l|4] . In magnetic films with pla- 
nar magnetization, the DW are generally large and are well 
characterized using classical tools such as magnetic force 
microscopy (MFM) [5 , scanning electron microscopy with 
spin polarization analysis (SEMPA) [6] or x-ray photoe- 
mission electron microscopy (X-PEEM) [7 . On the con- 
trary, in materials presenting a perpendicular magnetic 
anisotropy (PMA), the high magnetic anisotropy leads to 
narrow DW (width of few tens of nanometers or even less), 
which are difficult to measure and only the domains are 
observed [SHTO]. Based on micromagnetic arguments, in 
ultrathin films (thickness below the exchange length) the 
DW are generally thought to be Bloch walls [TOHlBj , which 
minimize the wall internal dipolar energy. However, this 
fact has never been explicitly proved by direct observa- 
tion. This lack of evidence is surprising as these materials 
are widely studied, since they present original phenomena 
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such as creep in small field [14] and since their DW are 
thought to be easily moved using spin polarized current. 

In this Letter, we present a study of DW in magnetic 
bilayers (two ferromagnetic layers separated by a non mag- 
netic spacer) with perpendicular magnetization. Such 
structures have been recently studied and the dynamics 
under field or current has been shown to be original due to 
the coupling between the two layers [15 . In order to gain 
precise information on their structure, we use a combina- 
tion of MFM and ballistic electron emission microscopy 
(BEEM), a recent technique that allows high spatial res- 
olution magnetization imaging. With the help of comple- 
mentary micromagnetic calculations, we show that, con- 
trary to intuition, the DW turn out to be Neel walls, due 
to the magneto-static coupling between the two ferromag- 
netic layers. Beyond this result, this study is also the first 
demonstration that BEEM can be applied to magnetic 
samples with perpendicular magnetization. 

Experimental details. — For the study of domain 
wall structure, high resolution magnetic imaging tech- 
niques have to be used. We have chosen to use the ballistic 
electron emission microscope (BEEM) [16], a STM based 
technique with a theoretical resolution of few nanome- 
ters JTl\. One advantage of this technique is that it is 
not restricted to samples grown and observed in a ultra- 
high vacuum (UHV) environment, a fact that, e.g. enables 
complementary magnetic force microscopy (MFM) mea- 
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surements. However, the BEEM technique also imposes 
some restrictions on the sample. In order to get enough 
signal, the magnetic layers have to be rather transparent 
to hot electrons in the 1-2 eV energy range. This forbids 
using Pt/Co multilayers, which is the archetypical exam- 
ple of multilayers with PMA, as the hot electron attenua- 
tion in Ft is very strong (more than 10 times higher than 
in Au [18 ). Instead, we use Au/Co multilayers, which 
provide similar magnetic anisotropics as Pt/Co multilay- 
ers but with a much better hot electron transparency. 

The samples have been grown using electron-beam evap- 
oration in UHV (base pressure 2 x 10~^^ mbar) on a 
hydrogenated n type Si(lll) substrate. The substrate 
has been cleaned and prepared using wet chemistry fol- 
lowing the procedure described in Ref. [19 , and dis- 
plays atomically flat terraces prior to metal deposition. 
The first 5 nm thick Au layer ensures the formation 
of a high quality Schottky diode (barrier height = 
0.8 eV, ideality factor < 1.1), which is a necessary con- 
dition for BEEM experiments. The following layers are 
Co(1.6 nm)/Au(5.0 nm)/Co(1.4 nm)/Au(5.0 nm) (the last 
layer is a capping layer which ensures the stability of the 
sample outside the UHV environment). The effective mag- 
netic anisotropy is expected to be rather low (but positive) 
due to the rather large Co thickness. The advantage is that 
the domain wall energy is expected to be small enough to 
allow for the nucleation of small domains, thus facilitating 
BEEM experiments. 

The BEEM experiments have been conducted using a 
modified Omicron-1 STM [20 at room temperature and 
in a UHV chamber (different from the evaporation cham- 
ber). A magnetic field up to 60 Oe and perpendicular to 
the surface can be applied from outside the UHV chamber. 
For a positive tunneling voltage Vt^ hot electrons are col- 
lected in the Si substrate and the transmission T defined 
as the hot electron current Ib to tunneling current It is 
measured. The MFM experiments have been conducted 
in ambient conditions on a Dimension 3000 Veeco AFM 
using home-made Co-Cr coated tips. 
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of the 

Co(1.6 nm)/Au(5 nm)/Co(1.4 nm) sample. (a) MFM 
phase image (inset topography) (b) BEEM (Ballistic current) 
image (inset topography) (imaging conditions: Vr = 1-5 V, 
It 20 nA) 



Results. — The sample have first been studied us- 
ing "classical" techniques, namely magneto-optical Kerr 
effect (MOKE) and MFM. The magnetic hysteresis loop 
measured in polar geometry shows that the sample indeed 
has a perpendicular magnetic anisotropy with a saturat- 
ing field of 500 Oe, but does not have a 100 % rema- 
nent magnetization [Mr /Ms ~ 20 %). This is an indica- 
tion that in zero magnetic field, spontaneous nucleation of 
magnetic domains occurs. This phenomenon, observed in 
many samples like garnets, is attributed to a low DW en- 
ergy, which allows minimizing the dipolar energy through 
the nucleation of small domains [21] . These domains have 
been revealed by MFM: the image (Fig. [TJi) shows two 
types of domains with opposite magnetization (perpendic- 
ular to the surface), each domain brand occupying 50% of 
the surface in zero magnetic field. The domain structure 
corresponds to segmented stripes with alternate magne- 



tization direction but with a rather disordered structure. 
The average width of the stripes is 350 ± 50 nm. As only 
two contrast levels are seen in the MFM image, the mag- 
netic state of the sample corresponds to either (i) one sat- 
urated layer and one layer with domains or (ii) a replica- 
tion of the domains in the two layers. BEEM experiments 
have then been performed as the two hypotheses should 
give very different results, since BEEM probes the local 
alignment between the magnetization of the two layers. 
In the first case, we expect images similar to MFM im- 
ages, whereas in the second case, we expect no contrast 
as the magnetization directions are always parallel to each 
other. 

BEEM imaging has been first performed in zero mag- 
netic field as shown in Fig. [TJd. The images are clearly 
different from MFM images but a clear contrast, uncor- 
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Fig. 2: Hot electron transmission (Ib/It) versus tunneling 
voltage, recorded on the bright areas and on the dark areas. 
The transmission raises above 0.8 V, which corresponds to the 
Schottky barrier height at the Au/Si interface. The contrast 
evaluated with the help of equation [1] is also presented. The 
tunneling current is It — 20 nA. 

related with the STM topography is observed: dark lines 
(low local ballistic electron transmission) with a width of 
few tens of nanometers are present. One important char- 
acteristic is that the lines never cross each other and thus 
delimit closed areas, which is typical of a border. As ex- 
plained before, such a result is unexpected. Assuming 
that the lines display the DW, these images correspond 
to a sample with a replication of the domains in the two 
layers: inside the domains the magnetization of the two 
layers are always parallel giving rise to a high transmission 
(bright areas). 

To quantify the transmission and the magnetic contrast 
of our sample we performed BEEM spectroscopy measure- 
ments and compared the curves registered over the dark 
features and over the bright areas (Fig. [2j). Below 0.8 eV, 
no ballistic electrons are collected as their energy is lower 
than the Schottky energy barrier (this threshold is the 
same in both regions, which proves that the contrast does 
not arise from defects at the interface). Above 0.8 eV, 
we can clearly see the difference of transmission between 
the bright and the dark areas. We define the magnetic 
contrast C (in percent) as 



C = 100 X 



Tdark 



(1) 



where Ti^right (resp. Tdark) is the transmission measured 
on a bright (resp. dark) area. This contrast is nearly 
energy independent, and is about 40 ± 10%. Using the 
hot electron attenuation length obtained in previous work 
for Co [ISlIlQ], the magnetic contrast for our sample is 
estimated to range between 40 % and 60 % at 1.5V. This 
good agreement proves that the magnetization of the two 
layers are indeed locally anti-parallel in the dark zones. 

In order to prove that these dark hues are DW, we have 
also performed BEEM imaging under magnetic field, ap- 
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Fig. 3: Evolution of the BEEM images with magnetic field 
(H). (a) STM topography, BEEM images at (b) H= Oe, (c) 
H= 60 Oe and (d) H= -60 Oe {Vt = 1.5 V, It = 20 nA, 
1x1 iiw?). The arrow in (b) indicates the profile shown in 
Fig.[6l 



plied perpendicularly to the sample surface. Fig. [3] com- 
pares ballistic current images recorded at the same posi- 
tion without magnetic field (b) and with a field of 60 Oe 
(c) and -60 Oe (d). Compared to H= Oe (Fig. [3b), 
the central domain is reduced for a positive magnetic field 
(Fig. [3b) while it is enlarged for a negative field (Fig. [3]i). 
Thus, it has been possible to move the dark lines without 
changing their width, which is typical of DW. Comparing 
BEEM images with the sample topography (Fig[3ti), it is 
possible to move one step further. Some defects (higher 
grains) can be observed in the topography and seem to act 
as pinning centers: the domain walls are "linking" some 
of these defects and only move in the defect free zones. 
This incidentally proves that the sample microstructure 
determines the coercivity. 

Discussion: origin of the contrast. — The origin of 
the anti-parallel alignment in the domain walls should now 
be explained. The interpretation is based on the dipolar 
coupling only, as the Au spacer thickness (5 nm) is suffi- 
cient to neglect any direct coupling [22 . We can rule out 
that domains are misaligned (Fig. [4^), which could give 
rise to an AP alignment near the domains border. Indeed, 
the width of the observed lines is very regular and field in- 
dependent (in the small field range < 100 Oe). Moreover, 
the stray field around the domain walls tends to stabilize 
the alignment of the domain walls as depicted in Fig. [4^. 
We thus rather claim that the dark lines correspond to 
domain walls having opposite chirality in the two layers, 
leading to the local AP alignment. In single ultrathin 
magnetic layers, DW are Bloch walls (helicoidal rotation 
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(a) Misaligned domains 
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(b) Bloch walls with opposite chiralities 
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Fig. 4: Possible configurations explaining the contrasts ob- 
served, (a) Misalignment of magnetic domains, (b) Bloch walls 
or (c) Neel walls with opposite chirality. The dashed lines 
sketch the flow of the stray field. 

of the spins in the DW plane) so that no magnetic charges 
are present, which minimizes the DW internal dipolar en- 
ergy. These walls may have two chiralities (left or right 
handed). However, as no specific stray field is associated 
with these walls, they cannot couple to each other (nei- 
ther ferro- nor antiferromagnetically ) . In our multilayer 
system, the probability to get opposite chiralities is then 
~ 50 %. Since the size of domains observed by BEEM 
is in agreement with MFM observations, we know that 
all the walls are characterized by an AP alignment of the 
magnetization of the two layers. 

Furthermore, it is easy to see that Bloch walls are unsta- 
ble in this multilayer sample. As shown in Fig. HJd, at the 
center of one wall the stray field of the domains in the other 
layer is aligned along the x direction (perpendicular to the 
wah plane). Since, at that position, the magnetization of 
Bloch walls is along the y direction, this field creates a 
torque, which destabilizes the magnetization toward the 
X direction. Due to this field, the walls may transform to 
Neel walls with a magnetization along the x direction, as 
shown in Fig. |4t. More interestingly, the torques devel- 
oped in each layer have an opposite direction, which leads 
to an opposite chirality in the two Neel walls. Moreover, 
as Neel walls are charged, an additional stray field, spe- 
cific to the walls, is induced, which even strengthens the 
AP alignment of the two layer magnetization in the walls. 
This effect is at the origin of the contrast observed in the 
BEEM images. Such configurations can be reproduced 
easily using micromagnetic calculations for magnetic lay- 
ers with a uniaxial perpendicular magnetic anisotropy 
with any value of parameters (i.e. thickness, anisotropy, 
magnetization...), provided that the magnetic anisotropy 



is larger than or equal to the shape anisotropy (|/io^|, 
with Ms the layer magnetization) and that the thickness 
is below the exchange length. 

We now try to get a quantitative analysis of our im- 
ages using numerical micromagnetic calculations. Us- 
ing the bulk Co values for exchange interaction (A = 
3 X 10~^^ J/m) and magnetization (Ms = 1-4 x 10^ A/m), 
the only unknown parameter is K which is difficult to 
measure, since our sample is not in a saturated state in 
zero magnetic field. We have performed 2D calculations 
[x^z plane), with infinite dimensions in the y direction 
(along the DW axis) for a system with stripes. In our 
in- house code [23], the sample periodicity in the x direc- 
tion is enforced. For a given the value of the period is 
then adjusted such as to minimize the energy density. The 
simulations reproduce well the model proposed in Fig. Hh 
with Neel walls with opposite chirality in the two layers 
(see Fig. [5^ and b). In order to get a quantitative com- 
parison, the goal is first to determine the value of K that 
reproduces the experimental periodicity P and compare 
the experimental and calculated DW profile. In a first ap- 
proximation, we consider two identical layers (same thick- 
ness, same magnetic anisotropy). In the range where K is 
of the same order of magnitude as the shape anisotropy, 
the period is a very good criterion, as P depends strongly 
on K as shown in Fig. [St. 

Experimentally, the periodicity (assumed to be at equi- 
librium) is 350 ± 50 nm so that, according to Fig. [St the 
quality factor is permitted to range (see Fig. [5t) between 
1.011 and 1.037 [K = 1.24 - 1.28 x 10^ J/m^). From 
the calculated magnetization profile, the BEEM signal can 
now be simulated. Supposing that the two spin conduc- 
tion channels in the spacer are independent, the BEEM 
signal is estimated to be [19] 

Tp + Tap , Tp - Tap ^ / n ^ / n /o\ 
Tbeem \ rni [x). m2 (x), (2) 

where mi {x) and 77?2 (x) are the local magnetization ori- 
entations in the top and bottom layers. Across the DW, 
the BEEM signal shows a deep minimum corresponding to 
the anti-parallel alignment of the magnetization. The DW 
width, defined as the full width at half minimum of the 
BEEM signal, is shown in Fig. [5b. For Q between 1.011 
and 1.037, we estimate a DW width of 40 ± 5 nm, close 
to the experimental value (43 ± 5 nm). In the experimen- 
tal BEEM images, we have recorded a profile in order to 
compare with the calculation as shown in Fig. [6] (note that 
several profiles on various images have been taken and they 
all show the same DW shapes). The very good agreement 
between the experimental and calculated profile strongly 
supports the magnetic structure model proposed here. 

In the sample, the actual Co layers thicknesses are 
slightly different (respectively 1.4 and 1.6 nm). To bet- 
ter describe the sample reality, we now consider two 
different layers in the simulations. Since the magnetic 
anisotropy has an interfacial origin, the total magnetic 
anisotropy of a film with a thickness t can be written as 
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Fig. 5: Micromagnetic calculation of a stripe network for two 
identical ferromagnetic layers (1.5 nm thick) separated by a non 
magnetic spacer (5 nm thick), (a) Schematic configuration of 
one stripe period (not to the scale) and (b) calculation result 
for K = 1.26 X 10^ J/m^ (c) Stripe periodicity P and DW 
width (see definition below) variation vs. quality factor Q = 
2K//ioM|. 
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Fig. 6: Comparison between the experimental profile (taken 
along the arrow in the BEEM image Fig [Sj)) and the calcula- 
tions. The profile has been taken through two DW and thus 
shows two minima. For the calculated profiles, the distance 
between the two DW has been adjusted manually after the 
calculation so that to correspond exactly to this experimental 



BEEM profile reproduces very well the experimental one, 
as shown in Fig. [6l 

Conclusion. — In conclusion, we have shown that in 
a magnetic multilayer with perpendicular magnetization, 
the usual Bloch walls turn out to be unstable. A new 
type of wall is found consisting of one Neel wall in each 
layer, with opposite chirality. The details of these walls 
have been revealed with a high spatial resolution using the 
ballistic electron emission microscope. The simulations 
reproduce very well the observed features and show that 
this type of wall occurs in such layers for any magnitude 
of magnetic anisotropics, consistent with the existence of 
a perpendicular magnetization. The consequences of such 
DW structure modification on their dynamics ought to be 
detectable. 



K = Ky + 2Ks/t, where Ky is the cobalt volume magne- 
tocrystalline anisotropy (0.4 MJ/m^ [24j) and Ks is the 
surface anisotropy. Thus, the thicker film is expected to 
have a smaller anisotropy than the thinner one. Consid- 
ering that the anisotropy found in the first calculations 
corresponds to the average of the anisotropics of the two 
films, we estimate Ks = 0.66 ±0.02 mJ/m^ in good agree- 
ment with values found from the literature [24] . The sim- 
ulation using Ks = 0.66 mJ/m^ and those different film 
thicknesses provides results almost identical to those pre- 
sented above, but with a larger DW width in the thicker 
layer. The equilibrium period is found to be 340 nm in 
good agreement with the experiments, and the simulated 
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